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BEHAVIOR OF THE IONIZATION CROSS SECTION IN HIGH 


ENERGY ELECTRON -HYDROGENIC ATOM COLLISION 
K . Qmidvar 

Laboratory for Theoretical Studies, Goddard 
Space Flight Center, Greenbelt, Maryland 

The ionization cross section per unit range of the 

energy e of the ejected electrons for collision of an 

incident charged particle of high energy E -with an 

atomic system is given by d Q/d e = E** 1 [A( e ) jgnE + B( e )], 

where A and B are functions of $, and the total cross 

section is given by Q = E 1 [ A inE + B ] . The coefficients 

t t 

A( e ) and B( e ) for the particular case of the atomic 
hydrogen are derived analytically as an expansion with 
respect to k = /g. The expansion converges for all 
values of k. By integrating A( e ) and B(e) with respect to 
e, A and B are obtained and compared with the values 

T T 

given by Bethe using a different method. The coefficients 
of the first two leading terms of the expansion of dQ/dg 
with respect to q or l/ e when e is s,mall or large Is found, 
lationship for the average energy of the ejected e. ctrons 
is given. Different features of the high energy imp ct 
ionization and certain connection between the ionization 
and the discrete level excitation is discussed. 


A re- 
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I. FORMULATION 

Bethe in a memorable articled has shown that for inelastic collision of 
a charged particle of high energy E with an atom the cross section dependes on E, 
and linearly on two parameters which are functions of the atomic structure and 
the energy transferred to the atom. For ionization of an hydrogenic atom of nuclear 
charge Z by an electron with non-relativistic high energy E, the cross section 
per unit range of the energy € of the ejected electron and the total cross section 
are given respectively by 

dQ/d e ~ (Z 4 k 1 2 )' 1 [A( e ) J enk 1 2 + B^Urra 2 , k x 2 = a 2 e/Z 2 , (l ) 

e ' 
max 

Q = J Q df dc r (z 4 k i 2 ) 1 [A t j^k , 2 + B T ] n a 2 ( 2) 

with a Q the Bohr radius and E and e in rydbergs. A( e ) and B( e ) are functions of 
e- A and B are dimensionless constants and independent of the nuclear charge, 

T T 

and are integrals of A( e ) and B( e ) with respect to e . e is the maximum value 

max 

of e . Eqs. (l) and (2) will be derived here. Bethe has derived (l) and (2), 
and using sum rules for the generalized oscillator strength he has evaluated A 

T 

and B for the hydrogenic atoms. Here A( e ) and B( e ) are derived analytically, 

T 

and by integrating them over e , values of A and B are rederived. 

T T 

With the system of an incident electron and an hydrogenic atom of nuclear 
charge Z let k* , k^, and k* represent the momenta of the incident, scattered, 
and ejected electrons, and K 1 = k* - k^. We take ti = m = 1 with m the mass of 
the electron, and throughout we express energy in rydbergs. Then E = k^ 2 and 
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e = k 1 2 . In this way E and e take the dimension of the inverse of the length 
squared, but numerically their values are given by rydbergs. In the Born 
approximation 


Q 


e 

J’max 

o 



(3) 


kj * 4 . 

d& - 8 tt r d|^_ df 

de a 2 E(Z 2 /a 2 + e) K' d e 


W 


where df/de is the generalized oscillator strength per unit range of e defined by 


df 

de 


Z 2 /a 2 + e 
' o 

2K 1 2 


J 


k' dk' 


<k* 


iK* 


X |ls>] i 


(5) 


with j ls> and | k' > the ground and the continuum wave functions of the atomic electron 


k 1 > is normalized such that asymptotically I k*> (2tt) 


ik 


We 


introduce the dimensionless vectors 

4, k 2 , k and K^by ^ = |Jk , kg = ^kg, k ' 1 = |nCj and K 1 = where n - ?./ a Q . 
Evaluation of df/de then leads to 1 * 2 


dQ _ 
de 


2 1 °( TT a 2 ) 

Z 4 E(l - e _2TT /k) 


k + kg 

1 

F(K)dK , 

k - kg 
1 


(6) 



_ k - 


f(k) 


e -g f -p/k j- K a + (•[_ + k g)] 

K[1 + 2(K 2 + k 2 ) + (K' 2 - k 2 ) 2 } 


T =? tan 


2k 


E 2 - k 2 + X 


(T) 


We seek an expansion of p(K) -which would converge for all values of the 
momentum transfer K. Similar expansion with useful application has been 
introduced by Lassettre 3 for excitation of the .non-hydrogenie atoms and 
deduction of the form of the differential cross section from experimental data. 
It is feasible that the present expansion will also find usefulness for 
ionization of the non-hydrogenic atoms. F(K) in (7 ) can be expanded in terms of 
K for small and large values of K, corresponding to the limits of the integral 
in (6) at high energy but the integrals of the two expansions may and, in fact, 
do differ by a constant. To overcome this difficulty we change the variable of 
integration from K to cp. This results in 


F ( K)dK H f 1 + ] dcp . 

2 ( 2 k) 5 • 2k cos cp - (X - k 2 ) sin cp 


( 8 ) 


To separate out the logarithmic term in the integral in (6) we introduce cp as 
the zeroth order term of the Taylor's expansion of cp with respect to K 2 . Then 
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cp i = 9 - «Po, 'Po = tan 


2k 


1 


k 2 


y 


(9) 


and by direct substitution 


JV(K)dK = 


-29o/k 


6(2k) 5 (1 + k 2 ) 


J e" 2Cpl / k sin 4 (fo + 9j ) [2k cot - 2 (l + 2k 2 )]dc? 1 . 


( 10 ) 


From (7) and ( 9 ) it is seen that for all values of K and k, cp x lies between 0 and 
- n . The transcendental functions in the integrand in (lO) can then be expanded in 
powers of cp x and the integration can be carried out term by term. With this procedure 
we obtain 


C9 

jF(K)dK = 


e" g3(p0 I - 4n(sin 9i ) + 1 + 2k 2 
6(1 + k 2 ) ^ (1 + k 2 ) 4 8 (2k) 5 


Re (3< 


-x 


4e 




4ic Pn -1 


+ e 


) 


(1 + k 2 ) 4 


r 


AjVi) 

m 


(-2 9 i ) 


m 


m =.l 


x 


m. 


], 


(n) 
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A m (9) 


- 1 ( l_Lkf)T (i _ ) f^ 211 ) R e (3 a ra - he 

“8 \ 2k ' L ^ ra L (m + 2n) t2nj X ^ 


XnTT 

» = 1 


^ 2190 ^ 


+ e 


4icp 


o m\ 
c ) 


1 * 2k g 
2k 


(ja 


m - 1 


. 2190,111-1 . 4i9o rn 
4e Id + e c 


■>i 


( 12 ) 


where p = l/k, a = p, b = p - i, c = p-2i, and B^ Is Bernoulli's number. 

Re stands for the real part of a complex quantity. 

In this way, through (6), dQ/dg is expressed as a converging sum over C p i 
for any value of the momentum transfer, and the problem of integration is 
reduced to a summation. For high energy incident electrons to first order, 

~ = (l t k 2 )/2k^, and k^ + k 2 = 2k ^ . Therefore the lower and the upper 

limits of the integral in (ll) are given at high energy by = k/2k^ 2 
and cp x = cp Q ^ respectively. Evaluating (ll) at the limits, combining it with 
(6), and comparing the result with (l) gives 


A(e) 


g y gB<p ° , 

3ti a (l - e“ 2nfi j (1 + k £ ) s 


(13) 


B(e) = A(e) [^A m (cpo) x | (l + k 2 )] . 

m = 1 


( 14 ) 
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Substitution of (l) into (3) and comparison -with (2) shows that 


CO 


A 

T 


AtcUs 


J A(- 2 )cOv' ? , 

O 


B( e )d, 


J B^Vdk 2 

o 


(15) 


where, because of the form of A( e ) and B( e ), the upper limits of the integrals 
have been extended to infinity without loss of accuracy 4 . Here, by definition, 
A(k 2 ) ■= A(e)d e /dk 2 - p 2 A( e ). Similarly B(k 2 ) = g 2 3( s ), and by (13) and (l4) 
A and B "become independent of the nuclear charge. "Numerical evaluation of 

T T 

the integrals in (15) yields 


A 1.15^65 

T 


B T = 5-02625. 


(15a) 


This completes the derivation of Eqs. (l) and (2). Using these values of A_^ 
and B in Eq. (2) we can write this equation alternatively as 

j 


4 TT a 2 M 2 

Q -2— #i4Ckf , M 2 = 0.28341 , C = Vo -047483 (2a) 

Z 4 k x 2 

Eq. (2a) is the form originally given by Bethe with M 2 = O.285 and 0 - Vo .048. 

A more accurate evaluation of M 2 and C using Bethe 1 s method has Been done by 
Inokuti 5 , and values given by (2a) are in agreement with these values to all 
places. Similarly Inokuti 6 has given the numerical values of A( e ) and B(e) in 
slightly different form, and (l5) and (l4) agree with these results. 
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II. RESULTS AND DISCUSSION •*, . 

It is convenient to define 'by dQo/dg the cross section -which is obtained 
by (V) by replacing the generalized oscillator strength by the optical 
oscillator strength, df^de, given by the first non vanishing term of the 
Taylor 1 s expansion of e ~ ~ in (5) . Then by straightforward integration. 


df o Ae = | ( Z 2 /a Q 2 + e) Jk' d£ | <k' |r|ls >| 2 




-2P9o 

e" 2TTP ) (1 + k 2 ) 4 ‘ 


( 16 ) 


Comparison of (l6) with (13) shows that 


A( e ) 


2 

3 




<k ! 


r 


ls> 


2 


(17) 


relating A(e) to the squared modulus of the dipole moment. Integrated over 
directions of the ejected electron. Substitution of (l6) in ( 4 ), replacement of the 
limits in the integral In (k) by their first order term, and an intergration 
leads with the use of (lj) to 


de 


2TrA(e) 

Z 2 E 



1 ■+ k 2 


(18) 


Comparison of this with (lj) and (ik) shows 


that (l) can be written 


dg ~ 1 dQ, o ttA ( g ) 
de 2 de + Z 2 E 


) A (9 ) 
^ m o 


m = 1 


(29 0 ) m 

in! 


X 


(19) 
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Then, contribution of the small scattering angle, or equivalently small 
momentum transfer by the incident electron, to the cross section represented 
by the logarithmic term in (l) is half the contribution of the dipole term 
integrated over all range of the momentum transfer. This provides a convenient 
way to avoid the introduction of a cut off in the integral in (k) for the 
distant collisions. 

The two terms in (19) correspond to the dipole and the non-dipole parts 
of the cross section, and they are comparable in magnitude throughout the 
non-relativistic region of the incident energy (Cf. Fig. 3). This is in 
contrast to the excitation of the low lying levels of the atoms where at 
energies several times the threshold energy the dipole part dominates the 
total cross section. This also indicates that while photo -excitation is closely 
related to the high energy electron impact excitation, such relationship does 
not exist between photo- ionization and the electron impact ionization. 

An interesting part of the present paper is to find the leading terms of 
the expansion of d Q/d e in terms of € or l/ e where e is small or large, respec- 
tively. Using the explicit form of A^Ccp) given by (12) we find after some 
algebraic manipulation that we have the following limiting cases when k is 
small or large 



1 ~ 2 r (4l e 4 + ll) + 2~7(l3e 4 


- H3/3 )k a . 


CO 
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22° (i * 3 ) 

2= k 3k= 


( 21 ) 


It also is necessary to expand A( s ) for small and large g. In doing so we can 
establish that for 0 k < 1 the factor (l - e 2TT / k ) -which appears in (15) 
is equal to -unity plus terms of the order less than k 17 . Then for the 
expansion considered here this factor can he put equal to unity and we obtain 


A ( e ) ** 
k -♦ 





( 22 ) 


k 


2 e 

3 n ^ 2 k 9 


C 1 



TT g " 5 

3 k 2 



Substitution of (20) through (23) into (l^) leads to 


B(e) ~ ^ (B - 3k 2 ), B = 4 • 62393, B 1 = 6.35310 

jj 2 O o 

k ~+ o 


(23) 


(24) 


p 2 k 4 


n - 


3k e 


(25) 



With these values of A( e ) and. B( e ) Eg., (l) can he written 


a 2 e 


dQ/de ~ (Z®^ 2 ) 1 [(A q - A x lc 2 ) m k 1 2 ' + B q - B x k 2 ] ua^, k 2 = ~ 

E » 1, e « 1 


A q -• 3.12587, A x =• ll.J1.a5 


cSA 

de 


1 « e « E 


4na 2 

o 

F.(a 2 ef 


0 -5f •) 


5*o * 


Eq. (.26) is a valid, expression when the incident and the scattered 
electrons have high, and the secondary electron has low, energies. Most 
contribution to the total cross section is due to this energy region of the 
secondary electron. Eq. (.26) is related to the excitation to a high principal 
quantum number n of the normal state of the hydrogen atom in the following 
way 7 . Let Q(n) be the excitation cross section. By writing e = - Z 2 n 2 a Q 2 , 
where e is the energy of the states with the principal quantum number n, 
and noticing that dQ/de is a continuous function of e when e becomes negative, 
through the relationship Q(n)dn = (dQ/dg)de we obtain up to n 4 terms: 

Q(n) ~ 2(n 3 Z 4 k 1 2 ) 1 [(A o + kjn z )*n k x 2 + (B q + B^/n 2 ) ]TTa Q 2 . 
n » 1 


Brown 8 has evaluated Q(n) for n » 1 by an independent method, and (28) 
agrees with his results. Similarly for high n and high incident energies 
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(28) is in agreement with the results of the numerical integration of Q(n) by 
the author 9 . 

Concerning Eq. (27) the leading term in this equation is the same as the 
elastic cross section between two electrons in a frame of reference where 
one of the electrons is initially at rest and its energy becomes q after the 
collision. The cross section in this case is the same classically and quantum 
mechanically, the classical cross section being given by the Rutherford formula 10 
and the quantum mechanical by the Moeller 1 s formula 11 . The reason for the agree- 
ment of the leading term with the elastic cross section is given in reference 10. 
Eq,. (27) in addition gives a second order term which is due to the binding of 
the atomic electron to the nucleus. This term agrees with the classical cal- 
culation of Vriens and Inokuti, where quantum mechanical momentum distribution 
is assumed for the bound electron (private communication). 

When e becomes comparable to E, the next higher term in E and also the 
contribution from the second Born approximation should be included in (27). 

In addition, the exchange effect for the identity of the two electrons should 
be taken into account. Because of the difficulty in calculating the second 
Born approximation this will not be dene here. Neglecting the second term in 
(27), the inclusion of the exchange effect in the leading term of (27) leads 
to the Mottls formula 10 


dQ ~ 
de 


IfTT 

a 2 E 
o 



1 

+ 



1 « Q, E 


(E - e ) 2 


e(E - e) 


(29) 
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Concerning the contribution of the exchange effect to the total cross 
section it should be pointed out that this is of the order of 1 /e, and 
therefore it becomes negligibly small at high enough incident energies* This 
can be seen by integrating (27) over the energy region -where q is comparable 
to E. 

. A quantity of interest is the average energy of the ejected electrons in 
an ionization process. This quantity -will be calculated below. By definition 


€ max e max 



(50) 


By writing e = |a 2 k 2 and by making use of (l) we find that 

k 2 - 1 k 2 - 1 

e ~ { J 0 [A(k 2 )jJnk 2 + B(k 2 )] dk 2 } { J Q k 2 [A(k 2 ) jjnk 2 + B(k 2 )]dk 2 } 

o 

(31) 


The integrals occuring in the denominator of ( 31 ) are given by (2) . When k is 
large A(k 2 ) and B(k 2 ) are proportional to k -9 and k" 4 respectively [Cf. Eqs. (23) 
and (25)]. The upper limit of the integral with respect to A(k 2 ) in the numerator 
of ( 31 ) can then be put equal to infinity. For the integral with respect to 
B(k 2 ) we have up to the first order terms 
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k 2 - 1 


k i 2 - 1 


J k 2 B(k 2 )dk 2 = J k 2 B(k 2 )dk 2 + k 2 dk 2 


a 

k 2 B(k 2 )dk 2 - ksx®. + - !)> (52) 

o 



where a is a value of k 2 at which B(k 2 ) reaches its asymptotic form. In this 
way we find that 


_ /Z , 2 C ^x 2 + D 

e = 

o A j^nk 2 B 

T ■* « 


( 33 ) 


with C and D dimensionless constants given by 



and A and B given by (ija) . 
t t 
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It is more appropriate to use an asymptotic form for B(k 2 ) in -which 
the exchange effect is taken into account . Comparison of the leading term 
of (25) with (29) shows that when exchange is included the asymptotic form 
of B(k 2 ) is given by 

+ — - — - - — 1 > 

L k 4 (k x 2 - k 2 ) 2 k 2 (k x 2 - k 2 ) J 

By making use of this form in (jl), and noticing that when exchange is included 

e = ~ ii 2 (k 2 - l), we find that (33) is still valid provided we replace D 
max <z_ 

in this equation by 

D ex = D + 4(1 - # 18 ) = - 6.it4l ( 35 ) 

It should be realized that although both the numerator and the denominator 
Of (35) are of the form j£hk 2 4. constant, the origin of the logarithmic terms 
is quite different. The term in the denominator comes entirely from small 
momentum transfer, that in the numerator partly from these and partly from 
large momentum transfer. This indicates a large spread in the energy of the 
ejected electrons. 

With or without exchange^ 7 as the incident energy increases reaches a 
limit given by (C/A ) (z/a ) 2 = ^-.063 Z 2 ryd. Then the maximum value of J 

T 0 

for an electron bounded to a central charge Z is proportional to its binding 
energy, although from (33) it is seen that the higher Z the higher the incident 
energy at which this maximum is reached. For many electron atoms it can be 
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said that the maximum value of the average energy of an ejected electron is 
roughly proportional to the square of the effective charge seen by the electron 
before its ejection; ; , 

We also notice that for atomic hydrogen the maximum energy loss* by the ' - 
incident electron cannot exceed 5..O63 ryd. Bethe 1 * 10 has shown that the 
energy loss by an electron is an inelastic collision with the “ hydrogen at cm 
reaches the value of 2 ryd at high energy, and that 28.3A$ of the collisions 
lead to ionization. Using these informations we find that at high incident 
energies the average energy loss of an incident electron in excitation of all 
the discrete states of the atomic hydrogen reaches the value of 0 . 7887 ryd. 

Figures 1 through A illustrate some of the results obtained. Fig. 1 is 
the well known graphical comparison of the Bethe approximation with the Born 
approximation, given here for the specific case of ionization. Z 2 EQ is plotted 
versus jgn(E/Z 2 ), where by Eq. (2) a straight line is obtained for the Bethe 
approximation. In the case of the atomic hydrogen the two approximations agree 
with each other , for energies larger than about 20 threshold units. For atoms 
with larger central charge the agreement is reached at still higher energies. 

Fig. 2 is a plot of the dimensionless coefficients A(k 2 ) and B(k 2 ) versus 
the dimensionless parameter k 2 = a 2 e/z 2 . The coefficient A(k 2 ) which is pro- 
portional to the dipole moment, Eq. (17)* is appreciable only for small values 
of k 2 , corresponding to the small momentum transfer by the incident electron 
arising from the distant collisions. For large k 2 it falls off as k 9 . The 
coefficient B(k 2 ) approaches its classical limit of A/k 4 given by the Rutherford 
formula. The figure shows clearly the transition between the quant al region, 
small e, and the classical region, large e* In the classical impulse approxi- 
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mation 12 16 the total ionization cross section at high incident energy is 
specified by ( 2 ) where A = 0 and B = 2 / 3 . 

T T 

In Fig, 3 the ratio of the logarithmic part of the cross section as 
given by (2) to the total cross section is plotted. It is evident that 
contribution of the dipole term given by the logarithmic term to the total 
cross section remains a fraction of the total cross section for the whole 
non-r elat ivis t i c range of the incident energies • 

In Fig. ^ a plot of the average energy of the ejected electrons versus 
the incident energies is given. The average energy reaches the limit of 
^.063 ryd as the incident energy tends to infinity, signifying the limit 
of the impulse received by th* atom. Relativistic corrections are neglected. 

To complete the discussion we list the coefficients A(n) and B.(n) originally 
given by Bethe for the electron impact excitation. Recalling the definition 
of Q(n) we can write in anology to (l) 

Q(n) - (Z 4 ^ 2 )" 1 [ACn )^ 2 + B(n)] n a 2 , = a 2 E/z 2 . ( 36 ) 

The dimensionless constant and charge independent coefficients A(n) and B(n) 
are given for n = 2 through n - 12 in Table I, 17 . For higher values of n Eq. ( 28 ) 


can be used. 
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Table I. Values of A(n) and B-fn) 


n 2 3 

4 

5 

6 

7 

8 

9 

10 

11 

n 3 A(n) 17-7577 9.6108 

7.9165 

7.2596 

6.9313 

6.7423 

6.6231 

6.5430 

6.4863 

6.4448 

n s B(n) 12 . 2 *'-' 10. 7054 10.0612 

9-7647 

9.6051 

9.5096 

9.4478 

9.4056 

9-3755 

9.3535 


12 

6.41J5 
9- 33ft 
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It is interesting to notice that the ration A/B for excitation and 
ionization is a monotonically decreasing function of the energy transferred 
to the atom, and at high energy transfer it approaches zero. The ratio, of 
the logarithmic term to the non-logarithmic term in the cross section is 
(A/B) j&hk^. For an incident energy of 1000 ryd and Z = 1, this ratio is 
10.022 for n - 2 excitation, and is I.558 for ionization. 

The foregoing results for electron impact ionization and excitation are 
also applicable to the proton impact ionization and excitation. This results 
from the fact that the cross section is proportional to the square of the 
charge of the projectile, and in the approximation stated following Eq. (12) it 
depends on the initial relative velocity of the projectile -target atom only, 
and is independent of the reduced mass of the system (Cf. Ref. 2). Thus all 
equations derived so far can be used for proton impact provided E is replaced by 
(m/M)Ep, where M and Ep are the mass and energy of the incident proton. In this 
way (36) agrees with the high energy analytical expression of Bates and Griffing 18 
up to, and excluding, the k 1 4 terms for n = 2 and n - 3 proton impact excitation; 
and (2) gives the same value at 3.5 MeV Incident protons as the graphical Bom 
calculation of Bates et al 18 , and Peach 19 , for proton impact ionization. 
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Figure Captions 


Figure 1. 


Figure 2. 
Figure 3 « 

Figure k. 


A comparison of the Be the 1 s approximation, Eq. (2), -with the . 
Born approximation* Q is the ground state ionization of an 
hydrogenic atom. 

A plot of the coefficients A(k 2 ), B(k 2 ), and 4/k 4 . 

is the logarithmic part of Eq. (2) and R is the ratio of this 
part to the total cross section. 

A plot of the average energy of the ejected electron versus the 
incident energy. Born approximation has been used for incident 
energies less than 25 ryd, and Eq. ( 53 ) for energies above this 


value . 
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